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We calculate the angular resolution of the planned LISA detector, a space-based laser interferometer 
for measuring low-frequency gravitational waves from galactic and extragalactic sources. LISA is not 
a pointed instrument; it is an all-sky monitor with a quadrupolar beam pattern. LISA will measure 
simultaneously both polarization components of incoming gravitational waves, so the data will consist 
of two time series. All physical properties of the source, including its position, must be extracted from 
these time series. LISA's angular resolution is therefore not a fixed quantity, but rather depends on the 
type of signal and on how much other information must be extracted. Information about the source 
position will be encoded in the measured signal in three ways: 1) through the relative amplitudes 
and phases of the two polarization components, 2) through the periodic Doppler shift imposed on the 
signal by the detector's motion around the Sun, and 3) through the further modulation of the signal 
caused by the detector's time-varying orientation. We derive the basic formulae required to calculate 
the LISA's angular resolution AQs for a given source. We then evaluate Afig for two sources of 
particular interest: monchromatic sources and mergers of supermassive black holes. For these two 
types of sources, we calculate (in the high signal-to-noise approximation) the full variance-covariance 
matrix, which gives the accuracy to which all source parameters can be measured. Since our results 
on LISA's angular resolution depend mainly on gross features of the detector geometry, orbit, and 
noise curve, we expect these results to be fairly insensitive to modest changes in detector design that 
may occur between now and launch. We also expect that our calculations could be easily modified 
to apply to a modified design. 
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I. INTRODUCTION 

This paper calculates the angular resolution of the 
planned LISA gravitational wave detector. LISA (short 
for Laser Interferometer Space Antenna) is in essence a 
space-based version of the ground-based interferometric 
detectors currently under construction: LIGO, VIRGO, 
etc. There are some major differences, however. The 
primary difference is that LISA will be sensitive to 
gravitational waves in a much lower frequency band: 
1CP 4 — 10 _1 Hz. (This low-frequency regime is unobserv- 
able by any proposed ground-based detectors, due to seis- 
mic noise. The ground-based interferometers will be sen- 
sitive in the range lO^lO 3 Hz.) The lO^-lO -1 Hz band 
contains many known gravitational waves sources that 
LISA is 'guaranteed' to see. These guaranteed sources 
comprise a wide variety of short-period binary star sys- 
tems, both galactic and extragalactic, including: close 
white dwarf binaries, interacting white dwarf binaries, 
unevolved binaries, W UMa binaries, neutron star bina- 
ries, etc. @-||. Indeed, our galaxy probably contains so 
many short-period, stellar-mass binaries that LISA will 
be unable to resolve them individually, and the resulting 
'confusion noise' will actually dominate over instrumental 
noise as the principal obstruction to finding other sources 
of gravitational waves in the datastream. Besides stellar- 
mass binaries, other possible LISA sources include: 1) 
a stochastic GW background generated in the early uni- 
verse, 2) the inspiral of compact, stellar-mass objects into 
supermassive black holes (SMBH's), and 3) the merger of 
two SMBH's [Q. The detection of any one of these would 
clearly be of immense interest. The events involving su- 
permassive black holes must surely occur in the universe, 
but the event rates are highly uncertain. 

Besides their different frequency bands, another im- 
portant difference between LISA and the ground-based 
interferometers concerns their means of identifying the 
angular position of the source on the sky. LISA is not 
a pointed instrument; it is an 'all-sky monitor' with a 
quadrupolar beam pattern. The ground-based detectors 
share this characteristic, but because there will be at least 
three ground-based detectors, and because they will be 
sensitive to gravitational radiation whose wavelength is 
much shorter than the distance between detectors, they 
will be able to determine the source position to within 
~ 1° by a standard time-of-flight method. This method 
is not available to LISA. Only one space-based detector 
is currently planned. Moreover the gravitational wave- 
length at the heart of the LISA band (~ 10~ 3 Hz) is of 
order 1 AU, so a second detector would have to be placed 
at least several AU away from the Earth for time-of-flight 
measurements to give useful constraints on source posi- 
tions. 

As we shall see, LISA can be thought of as two de- 
tectors, each measuring a different polarization of the 
gravitational wave. Thus the data consists of two time 
series. All information about the source position (as well 



as all other physical variables) must be extracted from 
these two time series. Angular position information is en- 
coded in the time series in the following ways. First, the 
relative amplitudes and phases of the two polarizations 
provide some position information. Second, most sources 
will be 'visible' to LISA for months or longer, so LISA's 
translational motion around the Sun imposes on the sig- 
nal a periodic Doppler shift, whose magnitude and phase 
depend on the angular position of the source. (In ex- 
actly the same way, radio astronomers take advantage of 
the Doppler shift caused by the Earth's rotation to deter- 
mine a pulsar's position to an accuracy much better than 
that implied by the beam width of the radio telescope.) 
Finally, as we shall describe, LISA's orientation rotates 
on a one-year timescale, which imposes a further source- 
position-dependent modulation on the measured signal. 
In Fourier space, the effect of the detector's changing 
orientation on a monochromatic signal of frequency /q 
is to spread the measured power over (roughly) a range 
fa ± 2/T, where T is one year. (The factor of 2 arises 
from the quadrupole beam pattern of the detector.) The 
effect of the periodic Doppler shift coming from the de- 
tector's center-of-mass motion is to spread the power over 
a range /o(l±u/c), where v/c ~ 10~ 4 . These two effects 
are therefore of roughly equal size at /o ~ 10 -3 Hz, which 
is near the center of the LISA band; rotational modula- 
tion is more significant at lower frequencies and Doppler 
modulation is more significant at higher frequencies. 

It is also worth emphasizing that much of the uncer- 
tainty in the position measurement arises from the fact 
that from this pair of time series one must try to ex- 
tract all the physical parameters of the binary: the or- 
bital plane of the binary, the masses of the bodies, etc. 
Errors in determining the source position are correlated 
with errors in these other parameters. The result is that 
LISA's angular resolution is significantly worse than one 
would suppose if one ignored these correlations. From 
this consideration, it should be clear that LISA's angu- 
lar resolution depends not just on the detector and the 
signal-to-noise, but on the type of source as well. 

In this paper we derive the formulae necessary for cal- 
culating LISA's angular resolution Afig for some given 
source, and we then evaluate Af2g for two sources of 
particular interest: monochromatic sources (e.g., short- 
period, white-dwarf binaries) and mergers of supermas- 
sive black holes. For these two types of sources, we per- 
force calculate how accurately LISA can measure all the 
other sources parameters, as well. A preliminary esti- 
mate of LISA's angular resolution has already been made 
by Peterseim et al. [@J|, but that estimate was only for 
monochromatic, high-frequency sources, and it assumed 
that the frequency, polarization, and amplitude of the 
signal were known a priori, so only the source position 
had to be extracted from the data. Also, for simplicity 
the estimate by Peterseim et al. H|| took into account 
the information from only a single polarization, and it ne- 
glected the information encoded via the rotation of the 
detector. In essence, our paper provides a much more 
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realistic calculation. Since LISA's angular resolution de- 
pends mainly on gross features of the detector orbit and 
noise curve, rather than exact details of the detector de- 
sign, we expect that our results will be fairly insensitive 
to contemplated design changes. We also expect that the 
calculation presented here could be very easily modified 
to apply to a different design. 

We now turn to our main motivation for considering 
SMBH mergers. The striking feature of these mergers is 
the huge amplitude of the emitted gravitational waves. 
(That is, huge compared to other GW sources!) LISA 
would be capable of detecting SMBH mergers at basically 
any reasonable redshift [z < 10, say) with signal-to-noise 
S/N ~ 10 3 , so long as the black holes are in the mass 
range 10 4 M© < M(l + z) < 1O 7 M . This mass range is 
set by the frequency band where LISA is sensitive. While 
the event rate for such mergers is highly uncertain (it 
could be several per year, or << 1/yr; see fl for a recent 
review), if SMBH mergers were discovered, they could 
provide a way of determining all the basic cosmological 
parameters - Hq, Oo, and Ao - to remarkably high accu- 
racy. The idea is that from the gravitational waveform 
one expects determine the luminosity distance Dl to the 
source to an accuracy of roughly (5/iV) -1 . (In fact, we 
find that LISA does roughly a factor 10 worse than this: 
typically ADl/Dl ~ 1%; see Table 2). As pointed out 
in the LISA Pre-Phase A Report 0] (hereinafter referred 
to as the LPPAR) , if the source position on the sky could 
be determined to sufficient accuracy that one could iden- 
tify the host galaxy or galaxy cluster, then presumably 
one could also determine the redshift optically. Clearly, a 
mere handful of such measurements would be sufficient to 
determine Ho, f^o, and Ao to roughly 1% accuracy. One of 
the key motivations for this paper is to see whether LISA 
has sufficient angular resolution to make such identifica- 
tions possible. In brief, we find that LISA will determine 
the SMBH location to no better than ~ 10~ 5 steradians 
(and typically to ~ 10 -4 steradians) which is not suffi- 
cient by itself to permit identification of the host galaxy. 
However this position measurement will be available days 
before the final merger; other telescopes (radio, optical, 
X-ray) should know when and where to look, so if the 
merger is accompanied by an electromagnetic outburst 
(e.g., from an accretion disk being carried along by one 
of the holes), the host galaxy might still be determined. 
Whether such an electromagnetic outburst can be ex- 
pected appears to be an interesting open problem jjj. 

The plan of this paper is as follows. Sees. II and III 
give brief summaries of relevant background information. 
Sec. II reviews the basic formulae of signal analysis and 
parameter estimation (mostly to establish notation and 
conventions), while Sec. Ill describes LISA's configura- 
tion and orbit, its response to gravitational waves, and its 
noise sources. In Sec. IV we derive LISA's angular resolu- 
tion for monochromatic sources, and in Sec. V we derive 
LISA's angular resolution for SMBH mergers. For both 
cases we calculate the Fisher matrix, which estimates how 
accurately the detector can extract all the physical pa- 



rameters of the system from the measured signal. Our 
conclusions are summarized in Sec. VI. 

We should state at the outset the principal limitations 
of this study. Firstly, and necessarily, the current de- 
tector design cannot be regarded as final, and the pub- 
lished noise curves-for both the instrumental noise and 
the confusion noise-can only be regarded as rough es- 
timates. Faute de mieux, we use the current design and 
the best estimates of the noise sources currently available. 
Relatedly, while it is a goal of the LISA design that the 
instrumental noise should be stationary and Gaussian [Q , 
probably we will not know how well this goal has been 
met until the the instrument is functioning. Faute de 
mieux, we assume the noise is stationary and Gaussian. 
Finally, to simplify things, in our treatment of SMBH 
mergers we assume that the two holes are in a circular 
orbit, and we assume that the plane of the binary orbit 
is fixed. That is, we ignore both orbital eccentricity and 
the precession of the orbital plane caused by the spins of 
the holes ||. We intend to take both these effects into 
account in a later paper ^(Tof. 

Throughout this paper, we assume that the detector 
arm length is much smaller than the gravitational wave- 
length, and we take the gravitational waveform h a b(t) to 
be in the (standard) de Donder gauge. Time is in units 
of seconds and frequency is units of Hz. 

II. REVIEW OF SIGNAL ANALYSIS 

This section briefly reviews the basic formulae of signal 
analysis, partly to fix notation. For a more complete 
discussion, see |ll]] or |12] ). 

The output of N detectors can be represented by the 
vector s a (t), a = 1,2, ...,N. It is often convenient to 
work with the Fourier transform of the signal; the con- 
vention we use is 

/oo 
e 2 * lft s a (t)dt, (2.1) 
-oo 

The output s a (t) is the sum of gravitational waves 
h a (t) plus instrumental noise n a (t). We assume that the 
noise is stationary and Gaussian. 'Stationarity' essen- 
tially means that the different Fourier components h a (/) 
of the noise are uncorrelated; thus we have 

(**«(/) Mf'T) = - f)S n (fU, (2.2) 

where '( )' denotes the 'expectation value' and S n (f) a p 
is referred to as the spectral density of the noise. When 
N=l (i.e., when there is just a single detector), we will 
dispense with Greek indices and just write s(f) and 
S n (f). 

'Gaussianity' means that each Fourier component has 
Gaussian probability distribution. Under the assump- 
tions of stationarity and Gaussianity, we obtain a natural 
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inner product on the vector space of signals. Given two 
signals g a {t) and k a (t), we define (g | k) by 



(g|k)=2/ [SM) 
Jo 



(2.3) 



In terms of this inner product, the probability for the 
noise to have some realization no is just 



p(n = n ) 



oc e 



-(n | n )/2 



(2.4) 



Thus, if the actual incident waveform is h, the proba- 
bility of measuring a signal s in the detector output is 
proportional to e -( s - h I s - h )/ 2 . Correspondingly given 
a measured signal s, the gravitational waveform h that 
"best fits" the data is the one that minimizes the quantity 
(s-h|s-h). 



It also follows from Eq. (2J3) that for any functions 
g a (t) and k a (t), the expectation value of (g|n)(k|n), for 
an ensemble of realizations of the detector noise n a (t), is 
just (g|k). Hence the signal-to-noise of the detection will 
be approximately given by 



* [h] = (h|h) 
N l J rms (h|n) 



(h|h) 



1/2 



(2.5) 



For a given incident gravitational wave, different real- 
izations of the noise will give rise to somewhat different 
best-fit parameters. However, for large S/N, the best-fit 
parameters will have a Gaussian distribution centered on 
the correct values. Specifically, let A* be the "true" val- 
ues of the physical parameters, and let A 4 + AA l be the 
best fit parameters in the presence of some realization of 
the noise. Then for large S/N, the parameter-estimation 
errors AA 1 have the Gaussian probability distribution 



p(AA 4 ) = Afe 



-I A.\ A A 



(2.6) 



Here is the so-called Fisher information matrix de- 
fined by 



r = 

1 y — 



ah 

9A 1 



d\i ) ' 



(2.7) 



and M = \/det(r/27r) is the appropriate normalization 
factor. For large S/N , the variance-covariance matrix is 
given by 



(AA'AA^) = (r 



0{S/N)- 



(2.8) 



III. LISA 

In this section we describe LISA's geometry and noise 
curve; these are the only aspects of the detector that are 
necessary for our analysis. For more details on how the 
detector works, we refer the reader to the LPPAR 



A. Detector Configuration and Orbit 

The geometry of the LISA mission, as currently en- 
visioned, is depicted in Figs. 1 and 2. The detector is 
a three-arm laser interferometer, with each arm being 
5 x 10 6 km long. It consists of six drag-free satellites, 
positioned so that two adjacent satellites sit at each ver- 
tex of an equilateral triangle. (One satellite at each ver- 
tex would suffice, but the current design calls for two, 
which provides some redundancy and greatly simplifies 
the pointing control.) The detector's center-of-mass fol- 
lows a circular, heliocentric trajectory, trailing 20° be- 
hind the Earth. The plane of the detector is tilted by 
60° with respect to the ecliptic; this angle allows the 
satellites to maintain the shape of an equilateral triangle 
throughout the orbit. We refer the reader to Fig. 3.17 
(p. 92) of the LPPAR for a useful picture of the orbital 
geometry, and to Faller et al. []l3| for a simple explana- 
tion of why the 60°-tilt allows the equilateral shape to be 
maintained. 

We label the arms 1,2,3, and call their lengths 
Li,L2,L^. Gravitational waves cause time-varying 
changes in arm lengths by amounts SLi(t), which are 
different in the three arms. The differences are measured 
interferometrically |14|. With 3 arms, there are two lin- 
early independent differences: SLi — 8L2 and 8L2 — 5L^. 
Therefore LISA will be able to measure simultaneously 
both polarizations of an incoming gravitational wave. 

We find it useful to introduce two Cartesian coor- 
dinate systems: 'unbarred' coordinates (x, y, z) tied to 
the detector and 'barred' coordinates {x,y, z) tied to 
the ecliptic. Unbarred and barred spherical polar an- 
gles (9, </>) and (6, <fi) are related in the usual way to 
the Cartesian coordinates: cos# = z/(x 2 + y 2 + z 2 ) 1 / 2 , 
cos# = z/ (x 2 +y 2 + Z 2 ) 1 / 2 , etc. We denote by x a the unit 
vector along the x-axis, and similarly for y a , z a ,x a ,y a ,z a . 
Here the superscript 'a' is an abstract index indicating 
that the object is a vector in 3-dimensional space. The 
unit vectors along the arms are called ij, 1%, if, respec- 
tively. The detector lies in the x — y plane, and the x — y 
coordinates rotate with detector. We assign, for all time, 
the following coordinate-directions to the three arms: 



I? 



cos7i x + sirrji y 



where 



7i = tt/12 + (i- l)7r/3 



(3.1) 



(3.2) 



as shown in Fig. 2. The 7r/12 term in Eq. ( |3.2| ) is in- 
cluded for later convenience. We choose the x — y plane 
to be the ecliptic (i.e., the plane of the Earth's motion 
around the Sun). The detector's center-of-mass follows 
the trajectory 



6{t) = tt/2, 0(i) = + 27rt/T 



(3.3) 



where T equals one year, and where 4>o is just a constant 
that specifies the detector's location at time t = 0. 
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The normal to the detector plane, z a , is at a constant 
60° angle to z a , and z a precesses around z a at a constant 
rate: 



z a -&[ cos<£(i) x a + sin <j){t)y a 



(3.4) 



Using z a lf = and Eqs. (|3.l| )-( |3~4"| ), we see that the direc- 
tions If can be written in terms of the barred coordinates 



as 



11 = COS CCi(t) 

+ sinai(t) 



cos 4>(t) y a - sin </>(t) x a 
f + i (cos 4>{t) x a + sin 4>{t) y a ) 



, (3-5) 



where the cti(t) increase linearly with time: 

ai(t) = 2nt/T - tt/12 - (i - 1)tt/3 + q , (3.6) 

where cyq is just a constant specifying the orientation of 
the arms at t = 0. 

In this paper we are primarily interested in LISA's an- 
gular resolution. The error box for the position measure- 
ment has solid angle Af2s given by 



2tt 



Afis A^s - (A/x s A^ s ) 



(3.7) 



where /is = cos 8 5 . The second term in brackets in 
Eq. 3.7 accounts for the fact that errors in fig and 4>s 



will in general be correlated, so that the error box on 
the sky is elliptical in general, not circular. Also note 
the overall factor of 2ir in our definition of Ailg. With 
this choice, the probability that the source lies outside an 
(appropriately shaped) error ellipse enclosing solid angle 
Afl is simply e -^V^ s ^ 




X 



FIG. 2. Illustrates the orientation of LISA's 3 arms in the 
x — y plane. 



B. Detector Response 

Because the LISA antenna has three arms, it produces 
basically the same information as a pair of two-arm de- 
tectors, and therefore is capable of simultaneously mea- 
suring both polarizations of the gravitational wave. To 
begin with, however, we shall consider only a single two- 
arm detector, formed by arms 1 and 2. The extension to 
two detectors will be straightforward. 



1, Single Detector 




FIG. 1. Shows the two coordinate systems used in our 
analysis: 'barred' coordinates tied to the ecliptic and 'un- 
barred' coordinates that are tied to the detector and rotate 
with it. 



We refer to the detector formed just by arms 1 and 2 
as 'detector I.' Detector I measures hj(t), given by 



hi(t) = [SL 1 (t)-SL 2 (t))/L 



\h ab {t)(l$l\-l%l§ 
2 \2 xx 2 yy) 



(3.8a) 

(3.8b) 
(3.8c) 



where we used Eq. ( 3A ) to go from ( 3 
is, for an arbitrary waveform h a b(t), the 



to ( [L8cl) . That 
relative arm 



length difference' hi(t) measured by a 60°-interferometer 
is always precisely -\/3/2 times as large as the 'relative 
arm length difference' measured by a 90°-interferometer 
(a result which is well-known among LISA cognoscenti). 
This result assumes the 60°-interferometer is 'placed 
symmetrically' inside the 90°-interferometer, as sho wn i n 
Fig. 2. That is the reason for our 7r/12 term in Eq. ([3.2]). 
This simple observation saves a lot of work. It means 
that all the formulae derived in the extensive literature 
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on 90°-interferometers (LIGO, VIRGO, etc.) can be car- 
ried over immediately to LISA: one just remembers to 
multiply the signal amplitude by v / 3/2. 

Consider a monochromatic, plane- fronted gravitational 
wave travelling in the — n a direction. The general such 
solution can be written as the sum of two orthogonal 
polarization states. Let p a and q a be axes orthogonal 
to n a , with q a = —e abc ni,p c . Define polarization basis 
tensors by 



H lb = PaPb - QaQb , H* b = p a q b + q a p b . (3.9) 

For a particular, unique choice of (p a ,q a ), called the 
waves' principal axes, there is precisely a ir/2 phase delay 
between the two polarizations: 

h ab (t) = A+H+ cos(2vr/i) + A x H* b sin(2vr/i) (3.10) 

Here A + and A x are constants, the amplitudes of the 
two polarization states, and we have omitted an arbitrary 
phase by our choice of the zero of time. Our convention 
is A + > \A X | > 0; A x > for right-hand polarized waves 
and A x < for left-hand polarization. 

The strain hj(t) that the waves produce in detector I 
depends on A + and A Xl the principal polarization axes, 
and the direction of propagation: 

hi(t) =^A + F+(6 s ,<fis,i>s)cos(2irft) 

+ &A x F?(e s ,<f>s,i>s)M^ft) , (3.11) 

where Ft and Ff are the "detector beam-pattern" co- 
efficients (l5|: 



\ (1 + cos 2 0s) cos 2cf> s cos 2ips 
— cos 6s sin 2(f> s sin 2ips , (3.12a) 



F+(0 8 ,<t>g,1> 8 ) 

F?{6 s ,<t>s^s) = 5(1 + cos 2 d s ) cos20 s sin2^s 



cos 9s sin 2</>s cos 2-05 



(3.12b) 



Here the subscript 'S' stands for 'source', (0s,4>s) give 
the source location in the 'unbarred', detector-based co- 
ordinate system, and ips is the 'polarization angle' of the 
wavefront, defined (up to an arbitary multiple of n) as 
follows: 



tamps = (z a q a )/(z h p b ) 



(3.13) 



In the case of interest here, the source polarization is 
assumed fixed, but the detector plane rotates through- 
out the observation. A very closely related problem 
was examined by Apostolatos et al. || (hereinafter re- 
ferred to as ACST), who investigated the case where 
the detector location and orientation are fixed (on the 
timescale of the observation) but where the orbital plane 
of the binary undergoes Lense-Thirring precession (due 
the post-Newtonian coupling between the bodies' spins 
and their orbital angular momentum), thereby modulat- 
ing the complex amplitude of the measured signal in a 
nearly equivalent way. We refer to ACST for derivations 



of many of the formulae quoted below. (There is one 
important difference between the 'rotating-source' case 
studied in ACST and the 'rotating-detector' case studied 
here: the 'Thomas precession' term identified in ACST 
is absent in the 'rotating-detector' case. Note also that 
some of the sign conventions in ACST arc different from 
those used here; in particular, A x is defined to be positive 
in ACST, but has no definite sign here.) 

To begin, we rewrite the signal (3.11) in the conven- 
tional amplitude-and-phase form. For a waveform whose 
amplitude and frequency are slowly changing functions 
of time, we can write 



hi(t) 



^Ai(t) cos| 



2Trf(t')dt' + (p pJ (t) + tp D (t)] 



(3.14) 



where fit) is the GW frequency that would be measured 
by a non-rotating detector positioned at the solar system 
barycenter, and Ai(t), tp p j{t), and are given by 



Mt)= (A%F+ 2 (t)+A%F? 2 (t)\ 
_ x /-A x F/(t) 



1/2 



r(t) = tan 



(3.15a) 
(3.15b) 



<p D (i) = 2nf (t)c" 1 R sm6 s cos ( ij>(t) - 4> s ) ■ (3.15c) 



where R = 1 AU. We refer to (fo(t) as the Doppler 
phase; it is just the difference betweem the phase of the 
wavefront at the detector and the phase of the wave- 
front at the barycenter. We ha ve neg lected second- 
order Doppler corrections to Eq. ( 3.15c ); this is justi- 
fied since such corrections are o f orde r (u/c) | ipo {t) \ ^ 
3 x 10~ 4 (//10~ 3 ) radians. In ( 3.15c ) we also neglect 
the small change in the source frequency / (as mea- 
sured at the barycenter) that occurs during the time de- 
lay -Rsin^s/c; the fractional correction to <fo(t) due to 
this effect is of order \f~ 1 {df /dt)R/c, which for a binary 
is - 0.04(4/^/M)(6Af/r) 5 / 2 (//10~ 3 ), where p and M are 
the reduced and total mass of the binary, respectively, 
and r is the orbital radius. 

Ai(t) is the waveform amplitude, and we refer to 
fp,i{t) as the waveform's polarization phase. The 
time-dependence of Aj(t) and <Pp,i(t) are determined 
by Fj~(i) and Ff(t), which in turn depend on 
0s(t),<f>s{t),a,ndif) S {t)- Using Eqs. (2.1)-(2.6) we find 

cos 6 s {t) = \ cos 9 S - ^ sin S s cos(0(f) - (3.16) 



?(*) 



tan 



7r/12 
V3 cos 6~s 



sin9 s cos(0(t) - (f> s ) 



2 sin 9s cos((f)(t) — 4>s) 



(3-17) 



Now consider the case where the monochromatic 
source is a circular, Newtonian binary. (This is quite 
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general: any monochromatic point source can be rep- 
resented as a circular, Newtonian binary.) The lowest- 
order, quadrupole approximation gives: 



e n b L c 
2M 1 M 2 
rD 



1 + {L a n a f 



rD 



(3.18a) 
(3.18b) 

(3.18c) 



where Mi and M2 are the two masses, r is their orbital 
separation, D is the distance between source and ob- 
server, and L a is the unit vector parallel to the binary's 
orbital angular momentum vector. The binary's circular 
orbit, when projected on the plane of the sky- i.e., pro- 
jected orthogonal to the waves' propagation direction- 
looks elliptical. The principal axis p a is just the major 
axis of this ellipse. (See_Fig. 21 in ACST.) We let L a 
point in the direction (9l,4 i l)- The angles 9s(t) and 
4>s(t) do not depend on the principal polari zatio n direc- 
tion p a , so they are already given by Eqs. (3.16)-( [3.17| ) 
above. Using Eqs. (3.13) and (3.18a), ips{t) is given by 



tan^ s (t) = ( L a z a - L a n a z b n b ) I ( e abc n a L b z c 



where 

L a z a = icos 6>l - ^sin 9 L cos (4>(t) - (f> 



(3.19) 
l ) (3.20) 



L a n a = cos 0£cos 9s + sm9p sin 6*5 cos — 4>s) 



(3.21) 



e a bcn a L b z c = ~ sin#L sin 9s sin (4> L - 4> s ) 

— ^cos(f>(t) ^cos^l sin 05 sin ^5 — cos #5 sin0£ sin^L^ 

— ^sin^(i) ^cos^5 sin#£ cosc/>l — cosOz sin^s cosfis^J . 

(3.22) 



To recapitulate, Eqs. (3.14)— (3.15) give hj(t) in terms 
of S , 4> s , A + , A^ , F + (t), and F x (t). A + and A x are 



given by Eq. ( 3.18 ), while the terms F + (t), an d fx ( t) are 
given in terms of 9s(t), <fts(~t), ^s(*) by E q. (ft. 12), an d 
finally s (t), 4> 8 (t\ are given by Eqs. (KlfQ ml), 

and (3.18)-(3.22). We note that an equivalent expression 
for Aj(t) was derived independently by Giampieri pq]. 



2. Two Detectors 

We have stated that with its three arms, LISA func- 
tions like a pair of two-arm detectors, outputting two lin- 



early independent signals: si(t) = 6Li(t) — 8L 2 {t) I / L 



and sip = \^5L 2 (t) — SL 3 (f)j / L. We now extend our 

above analysis to include both detectors. We continue to 
assume the noise is stationary and Gaussian; nevertheless 
the noise in different arms will generally be correlated. 
Let rii(t) = SLi(t)/L be the noise in the i th arm; then 



(n i (m j (f')) = C ij (f)S(f-f) 



(3.23) 



where Cjj (/) ^ in general. Clearly then the noise out- 
puts of detectors / and //' will be correlated too. How- 
ever it is also clear that we can always find some non- 
trivial s/7 which is a (frequency-dependent) linear com- 
bination si and sn', and which is orthognal to si, in the 
sense that the noise in sn is uncorrelated with the noise 
in si. Just set 



hxV) = hi>(J)-Hf) 



C 12 (f)-C 13 (f)-C 22 (f) + C 23 (f) 



Cn(/)-2C 12 (/) 



and then we have 



c 22 (.f) 

(3.24) 
(3.25) 



We find thinking in terms of such orthogonal detectors 
to be very convenient for calculations. 

Unfortunately, there do not yet exist estimates for how 
the noise in LISA's three arms will be correlated. In this 
paper we will make the simplifying assumption that the 
noise is 'totally symmetric' among the three arms, in the 
following sense: 



C 12 (f) = C 23 (f) = C 31 (f), 
C n (f) = C 22 (f) = C 33 (f) . 



(3.26a) 
(3.26b) 



It is not too unreasonable to suppose that the instru- 
mental noise will approximately totally symmetric, since 
the individual satellites will all be virtually identical. 
Also, one can easily show that the confusion noise due 
to an isotropic background of gravitational wave sources 
must be totally symmetric. The reason is that for an 
isotropic background Cij(f)/ S n {f) can only be a func- 
tion of (e a bZf/^) 2 , where e a b is the Euclidean 3-metric. 
The totally symmetric condition then follows from the 
facts that LISA's arms are all the same length and the 
angle between any two different arms is 7r/3. (Of course, 
our galaxy is not spherically symmetric, and so the con- 
fusion noise from galactic binaries cannot be expected to 
satisfy the condition of total symmetry.) 

For totally symmetric noise, the linear combination of 
§i(f) and sn'(f) that is orthogonal to s/(/) is actually 
frequency- independent. That combination is 



sii 



(t) = 3~ 1/2 (sLi(t) + SL 2 (t) - 2 SL 3 {t) ^ /L 
= 3- 1/2 (si(t) + 2sn,(t^j (3.27) 
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which implies 

(n a (f)\Mf))=S a pS(f-f)S n (f), (3-28) 

where a,/3 take on values / or II, and S n (f) is the spec- 
tral density for detector I. 

In terms of the 'unbarred' coordinates introduced in 
Sec. III. A, it is easy to show using Eq. (3.1) that 



hn(t) 



2 "-xy 



2 y x J 



(3.29) 



That is, just as hi(t) is equivalent to the response of of 
a 90°-interferometer (modulo the overall factor V3/2), 
so hn(t) is equivalent to the response of another 90°- 
interferometer, rotated by 7r/4 radians with respect to 
the first one. (This result waspreviously derived by P. 
Bender, in unpublished work p7| .) It is therefore trivial 
to write down the beam coefficients for detector II: 



F+(6 s ,<j> s ^s) = F+(6 s ,<t>s 
F*(9 s ,<l> 3 ,il>s) = F*(9s,<l>s 



tt/4, ip s ) 



(3.30a) 
(3.30b) 



Finally, in complete analogy with Eqs. ( 3.14 ) and ( [3.15 ), 
we can write hu(t) as 

hn(t) = ^ A H (t) cos[ J 2nf{t')dt' + <p p ,n(t) + <p D (t)] . 

(3.31) 



where 

A H (t) = 



A 2 F + ' 



(t)+AlF*\t) 



1/2 



<P P ,ii(t) = tan 



-A.F*(t) 
A + F+(t) 



(3.32a) 
(3.32b) 



and where </>rj(i) is still given by Eq. (3.15c). 

Figs. 3 and 4 show an example of the modulation pat- 
terns due to detector rotation, for one representative 
choice of parameter values. 



FIG. 3. 
The amplitudes Aj(t) and Au(t) during a one- year obser- 
vation, for the following choice of initial detector position 
and orientation and of the source's position and polarization: 
0o = 0, Q = 0, fi s = 0.3, cf> s = 5.0, fiL = -0.2, 4> L = 4.0. The 
overall scale is arbitrary. 



4>p 




t (10 sec) 

FIG. 4. The polarization phases if> p j(t) and tp P ,ii(t) dur- 
ing a one-year observation, for the same parameter values as 
in Fig. 4. 



C. Noise Spectrum 



0.80 



A , , A„ 



0.60 




We will assume the noise spectra for detectors I and II 
are the same; both are given by S n (f), which we repre- 
sent as the sum of instrumental noise S n ^ n (f) and con- 
fusion noise tSV^cc^/)' 



Snij) — Sn,in(f) + <Sn,co(/) 



(3.33) 



These two contributions are shown in Fig. 5. We now 
consider them in turn. 



1. Instrument Noise 

It is a LISA design goal that the instrumental noise be 
stationary and Gaussian; our analysis will assume that 
goal has been met. The following is the current best 
estimate of the spectral density of the instrumental noise 
S n ,in(f) for detector I f§. 



SnMf) = 5-049 x 10 s 



a 2 (/) + /3 2 (/)+7 2 



(3.34) 



0.0 1.0 2.0 3.0 4.0 



t (10 7 sec) 



8 



where 
«(/) 



10 
10 



-22.79 (//1Q - 
-24.54 ^/ 1Q - 



3 )" 7/3 i 

5 ), 7 = 10 



-23.04 



(3.35) 



Here a(f) is mainly due to temperature fluctuations, 
(3(f) reflects the loss in sensitivity when the gravitational 
wavelength becomes comparable to or shorter than than 
the detector arm length, and 7 (a constant) is mai nly du e 
t o ph oton shot noise. We will assume that Eqs. ( |3.34| )- 
( |3.35 ) above give the noise spectral density for detector 
II as well. 



2. Confusion Noise 

It seems very likely that gravitational waves from 
short-period, stellar-mass binaries will actually be more 
important than the instrumental noise in 'drowning out' 
the signal from other types of sources. The issue of why 
stellar-mass binaries should be regarded as effectively a 
noise source is a subtle but important one. We confine 
ourselves to a few remarks on this subject, and refer to 
the reader to Bender and Hils [jl9| and Hils, Bender, and 
Webbink || for further details. 

The first remark is that the orbits of these stellar-mass 
binaries can be treated as Newtonian, and the radiation 
computed accurately from the quadrupole approxima- 
tion. The orbit of a Newtonian binary is periodic, so 
in Fourier space its gravitational radiation is composed 
of discrete lines at / = 2/P, 4/P, 6/P, , etc., where 
P is the orbital period. (The sequence is multiples of 
2/P instead of 1/P because the radiation is quadrupo- 
lar.) It seems reasonable to assume that most of the 
short-period, stellar-mass binaries have small eccentric- 
ity: e < 0.2. In this case, > 60% of the power comes 
out in the lowest-frequency line, f — 2/P [[|. Cur- 
rent estimates of the confusion noise therefore neglect 
the higher-frequency lines. 

Second, for an observation time of To ~ 1 yr, the dis- 
crete Fourier transform sorts monochromatic signals into 
frequency bins of width A/ = 1/T ~ 3 X 10~ 8 Hz. The 
typical timescale on which these binaries evolve is > 10 7 
yrs; thus in one year's observation, a binary's emitted 
GW frequency changes by < f/10 7 = 10~ 10 (//10- 3 ) 
Hz — i.e., much less than the width of one bin! Thus, ig- 
noring for the moment the motion of the detector, each 
binary remains in the same bin throughout the observa- 
tion. The lower half of the LISA band, 10~ 4 - 10~ 3 Hz, 
contains roughly 3 x 10 4 frequency bins, while our galaxy 
contains ~ 3 x 10 7 close white dwarf binaries (CWDB's) 
with GW frequencies in this range, so roughly 10 3 per 
bin. Thus the problem of 'fitting' for all the binaries, in 
order to then 'take them out' of the data, is extremely un- 
derdetermined: there are at least 10 3 times as many free 
parameters as there are data points. In fact, to model 
this signal accurately, one needs not only the frequency 
of each source but its location and orientation, since the 



motion of the detector 'smears out' the the signal over 
a frequency range ~ 2/year, in a manner that depends 
on these additional variables. So the motion of the de- 
tector only aggravates the problem of 'fitting out' the 
stellar-mass binaries. One might even argue that instru- 
mental noise is in principle no different from binary con- 
fusion noise; instrumental noise always arises from some 
deterministic physical processes that one could also in 
principle model and then attempt to remove from the 
data by some fitting procedure (or by monitoring these 
processes directly), but in practice one reaches a point 
where there are simply too many variables-too many free 
parameters-to obtain a fit that has predictive power. 

The number of galactic binaries per bin decreases with 
increasing frequency. Somewhere in the range 10 -3 Hz 
to 4 x 10~ 3 Hz, there is a transition from having many 
galactic binaries per bin to having fewer than one, on 
average. At frequencies above this transition, most of the 
information about some broad-spectrum source (such as 
a SMBH merger) will come from the bins that do not 
contain galactic binaries. At these frequencies, then, the 
binary confusion noise is dominated by the extragalactic 
contribution. 

The following represents the current best estimate by 
Bender and Hils [jl9| for the level of the binary confusion 
noise: 



Sn,co{f) 



10 


-42.685 f 


-1.9 


/<10" 


-3.15 


10 


-60.325 f 


-7.5 


10 -3.15 


< /'< io- 2 - 75 


10 


-46.85 j- 


2.6 


10 -2.75 


</ 



(3.36) 



where / is measured in Hz. This estimate assumes a 
space density of CWDB's which is 10% of the theoretical 
value predicted by Webbink ^0|, and does not yet take 
into account the contribution from helium cataclysmic 
variables (which are likely to be important in the range 
1 to 3 mHz). 
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-34.0 
log S n (f) 
-36.0 

-38.0 

-40.0 

-42.0 





-44.0 



-4.0 -3.0 -2.0 -1.0 

log f (Hz) 

FIG. 5. Shows the spectral density of both the instru- 
mental and confusion noise. log 10 S n ,m(/) and log 10 S , „ lCO (/) 
are plotted versus log 10 f. The total noise is given by 
(/) (/) + 



IV. MONOCHROMATIC SOURCES 

Any monochromatic point source of gravitational 
waves can be thought of as a circular-orbit binary: they 
are in 1-1 correspondence. Therefore, in the rest of this 
section we shall always speak of the angular resolution of 
LISA for circular-orbit binaries, but the results would ap- 
ply to any monochromatic source. An arbitrary circular- 
orbit binary is described by seven free parameters: the 
frequency fo (as measured by an observer at the solar 
system barycenter); the angles 6$, 4>s, 0l, and 4>l] an 
overall amplitude proportional to A = MiMz/rD; and a 
trivial overall phase fo related to the choice of t = 0. 



A. The Measured Signal 

For a circular, Newtonian binary, the waveform h a (t) 
(a — I, II) can be written as 



h a (t) = A a (t) cos Xa(t) 



(4.1) 



where 



(t) + <p D {t). (4.2) 

where tpo is just a constant of integ ration , and A a ( t), 
<Pp,a(t)> and ifiD{t) are given by Eqs. ( 3.1 5| ) and ( 3.32| ). 



The calculation of the Fischer matrix is simplified by 
the following trivial observation. Although the measured 
frequency is not exactly constant, due to the motion of 
the detector, it is very nearly the constant fo- Therefore 
we can take the factor 1/S n (f) out of the integral in 
Eq. (2J3), and write 



S„(/ ) 



Sn(/o) 



E 

a=I,II 

E 

i=i ji 



d t h* a (f)d 3 h a (f) df (4.3a) 
dih a (t) djh a (t)dt . (4.3b) 



where we used Parseval's theorem to go from (4.3a 



to 



(4.3b). Then using the fact that / >> A^dA/dt, we 
can approximate (4.3b) by 



(dihidjh) = sw(/ )- 1 £ / 

a=IJI J- 



A 2 a (t)diXa(t)d jXa (t) 



diA a (t)djA a {t) 
dt. (4.4) 



Thus to evaluate the Fisher matrix (4.4), we need the 
derivatives of A a (t) and Xa (t) with respect to the seven 
physical parameters In A, tpo, fo, 4>s, @L, 4>l- Clearly 
one might strai ghtfor wardly use the chain rule with 
Eqs. ( 3.15 ) and ( 3.32 ) to determine the partial deriva- 
tive^ of A a (t) and Xa (t) with respect to the four angles 
§s, 4>s, ^i,and <t>L, though the final expressions would be 
cumbersome. In our calculation, we preferred simply to 
take these derivatives numerically. The remaining partial 
derivatives are: 



dA a (t) 
d\nA 

din A 



A a (t). 



dA a (t) 
cVo 



= 0, 



dAJfi 
dfo 







d(fio ' dfo 



2-Kt 



(4.5a) 
(4.5b) 



B. Results 

For monochromatic sources, the detailed shape of the 
noise curve has no bearing on the Fischer matrix; all 
that matters is S n (fo), where fo is the source frequency. 
And S n (fo) is inversely proportional to the signal-to- 
noise of the detection, so we can eliminate S n (fo) from 
the problem simply by normalizing the results to some 
fixed, fiducial S/N. In Table I, we normalize our results 
to S/N — 10, where S/N is the total signal-to-noise accu- 
mulated by both detectors I and II. The advantage of this 
way of representing our results is that Table I remains 
valid for any noise spectral density. (However the results 
do depend on our assumption that the noise is totally 
symmetric among the three arms.) Table I lists LISA's 
angular resolution AQ$ f° r a one-year observation, for 
a range of source frequencies fo and for representative 
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choices of angles 9s, 4>s, 4>l- For these case, we also 
list Si/N and AO^/, the signal-to-noise and angular res- 
olution of detector I taken alone. (One expects Si/N to 
be roughly 2~ 1 / 2 S/N — 7.07, but the exact value clearly 
depends on the various angles specified.) The sizes of the 
position error ellipses Af2s and A fig / simply scale like 
(S/N)-\ 

Since everything in the problem is periodic with period 
one year, one obtains exactly the same results for T years 
of observation, when T is an integer. (That is, one ob- 
tains the same results after normalizing to S/N — 10; if 
instead one normalizes to sources at some fixed distance, 
then S/N scales like VT and Af2 s scales like T _1 .) This 
scaling will hold approximately, but not exactly, when T 
is some non-integer greater than 1. LISA's angular res- 
olution would certainly be much worse for observation 
times significantly less than one year. 

The results in Table I are easily summarized. LISA's 
angular resolution for monochromatic sources is roughly 
in the range 10 -3 to 10 _1 steradians (equivalcntly, 3 to 
300 square degress) for source frequencies in the range 
10 -4 < /o < 10~ 2 and S/N = 10. Having data from 
both detectors I and II provides hardly any improvement 
in angular resolution, apart from the trivial improvement 
due to the increased S/N. Presumably this is because in 
one year's time, LISA's changing orientation allows de- 
tector I by itself to measure both polarizations of the 
incoming wave fairly accurately. LISA's angular resolu- 
tion is roughly a hundred times better at 10 -2 Hz than 
at 10 -4 Hz (for fixed S/N). Clearly this is because the 
Doppler modulation is a much bigger effect in the higher- 
frequency sources. 

We note that one application of the results in this sec- 
tion is to compact stellar mass objects (like NS's) in 
orbit around SMBH's. In that case, orbital evolution 
is sufficiently slow on the timescale of the observation 
that the signal is effectively a sum of 'lines' at frequen- 
cies n/P, where P is the orbital period. To the extent 
that the strongest line dominates, the source is therefore 
monochromatic . 



V. SUPERMASSIVE BLACK HOLE MERGERS 

In this section we consider the information LISA could 
extract from the collision of two supermassive black holes. 
Note that since high signal-to-noise ratios are expected 
for this case, we expect the Fischer matrix approach to 
error estimation to be quite accurate. 

SMBH collisons are related by simple re-scaling to the 
mergers of stellar-mass black holes; many of the issues 
that arise in GW data analysis are identical for the two 
cases, so we take take advantage of the extensive liter- 
ature on the latter. Our exposition in this section will 
therefore be highly abbreviated; we refer to jllj for more 
extensive discussion and derivations. 

Unlike our treatment of monochromatic sources in 



Sec. IV, our treatment here will be based on two simpli- 
fying assumptions that probably are not justified physi- 
cally. Most importantly, we neglect the spin-induced pre- 
cession of the binary's orbital plane, and its correspond- 
ing modulation of the waveform amplitude and phase. 
Given that the rotation of the detector and the rotation 
of the binary's orbital plane modulate the signal in very 
similar ways, it may in practice prove difficult to dis- 
entangle these two effects in the data analysis, and the 
accuracy of position measurement would be correspond- 
ingly degraded. The second complication we omit is the 
possible eccentricity of the orbit. We plan to explore 
these two complications in later work. 

One might expect such collisions to occur at significant 
redshift (z > 1) so in this section we carry along the fac- 
tors of (1 + z) that we have ignored in earlier parts of the 
paper (see Ref. 0|). These collisions would be visible 
to LISA for a significant fraction of year. Assuming cir- 
cular orbits, the time interval from the instant t(f) that 
quadrupole frequency sweeps past / until the instant t c 
when the two bodies collide and merge is (to lowest or- 
der): 



t c -t(f) = 5(87rf)- 8 / 3 [M(l + z) 



-5/3 



= 3.003 x 10 6 s (//10 



_4^-8/3/ M (1 + 2) 
10 6 M o 



-5/3 



(5- 



where M = M 3/5 M 3/5 /{M 1 + M 2 ) 1/5 , and / and t are 
the frequency and time measured by an observer at the 
Earth. 



A. The Measured Signal 



We define a signal H(t) by 



H(t) 



2 Mi M 2 (l+z) 
D L r(t) 



/(tO df 



(5.2) 



where f(t) is the (redshifted) frequency a detector would 
measure if it were nonrotating and its p osit ion were fixed 
at the Solar system barycenter. In Eq. (5.2), Mi, Mq. and 
r are the unrcdshifted masses and orbital separation that 
would be measured by an observer near the source, and 
Dl is the 'luminosity distance' to the source |2l]]. We 
think of H (t) as a 'carrier' signal, which is modulated by 
the motion of the detector. The measured signal h a (t) is 
given by 

h a (t) = A a (t) cos (^J f(t') dt' + (t) + <PD(t) 

(5.3) 



where A a (t), defined by 



Aa(t) 



D L r(t) a ,,\ 
2 Mi A/ 2 (l+*) ' 



(5.4) 
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basically encodes the amplitude modulation. How is 
h a (f) related to H(f), the Fourier transform of the 
carrier? Since A a (t), tp Pia (t), and <fD{t) all vary on 
timescales of ~ lyr >> 1//, we can approximate h a (f) 
using the stationary phase approximation. This gives g] 



t{f)^t c -h{^fr^[M{\ + z)\ 



-5/3 



M/)=A Q (i)e-H<^ 



,(*) + ¥>£,(*) 



h(f) (/ > 0) (5.5) 



where i = t(f) is the instant at which the GW frequency 
sweeps through the value /. So given H(f), we have a 
simple (approximate) expression for h a (f). 

In this p aper, we use a 'model' of H(f) that was de- 
veloped in |22,0]. This model is based on the so-called 
'restricted PN approximation'. In brief, the model in- 
cludes post-Newtonian corrections to the phase of the 
waveform through 0(v/c) 3 , but it takes the waveform to 
be quadrupolar, with an overall amplitude given by the 
lowest-order approximation. A more detailed explana- 
tion of this model is given in |llj] . Specifically, our model 
H(f) is 



A f-r/6^ 0</</ c . 

f > j 'cut-off 



H{f) = l~ J ° »^J^tcut-off ( 5 _g) 



where 



A = {h/mfl^-V^Dl 1 [M(l + z)] 5/6 cA) 
*(/) - 2vr/i c - c»c - vr/4 + |(8tt.M(1 + z)f)~ 



(5.7) 



-5/3 



20 /743 Uu\ lin „ N ,„ 

1 + yU6 + i^r +(4/? - 167r)x 



Here M = M x + M 2 , (J, = MiM 2 /M, M = fj?/ 5 M 2 / 5 , 
and the PN expansion parameter x(f) is defined by 



x(f) = nM(l + z)f 



2/3 



(5.9) 



The parameters t c and <p c are just constants of integra- 
tion. The term (3 is a P 15 N spin-orbit coupling term de- 
fined in ((3 is only approximately conserved by the 
P 1 ' 5 N equations of motion, but in our model we treat 
it as a constant). In this paper, we 'cut-off' the signal 
(somewhat arbitarily) at f cu t-off = (3 3 / 2 7rM(l + z))^ 1 , 
corresponding to r/M — 3. In fact, the value of this 
cut-off substantially affects the calculated signal-to-noise 
S/N in cases where M(l+z) is greater than ~ 3x 1O 6 M0, 
because in these cases most of the S/N accumulates 
at the very end of inspiral. That is because LISA's 
noise curve S n (f) falls very rapidly with increasing / for 
/ < 3 x 10 -3 Hz. However we will find that the predicted 
angular resolution of the detection, for a source at fixed 
distance, is r athe r insensitive to fcut-off- 
From Eq. (5.5) we therefore have 



hj(f) = A(t)e t ( x{f) - p '' {t) - pD{t) )h(f) (f > 0) (5.10) 
where t = t(f) is given through 0(v/c) 3 by jll| 



4 /743 UA 327T 3/ 2 , 2) 
3 l 336 + 4A/J X 5 + 1 ' 



(5.11) 



Note that our model of the signal is the just inspi- 
ral waveform: it does not include the final merger and 
ringdown. Partly this is necessitated by our current ig- 
norance about the final merger. However it also seems to 
us that this neglect is justified by the particular prob- 
lem we are trying to solve; that is, it seems unlikely 
that signal from the final burst will significantly improve 
LISA's angular resolution, even if it dominates the signal- 
to-noise. The reason is that the final burst is only a 
few seconds long, and obviously LISA's orientation and 
velocity hardly change in that interval, so there's mod- 
ulational encoding of the source position. (Of course, 
the final burst could give extra information about the 
chirp mass M, which is correlated with the position un- 
knowns, but M. is already very well determined by the 
lower- frequency data.) 

Note that h a (f) depends on 9 physical parameters: 
M, fi, (3, 4> c , t c , laD L , 6s, 4>s, 4>l- The next step 
is to evaluate the Fisher matrix ( |2.8| ). As in Sec. IV, 
we numerically evaluate the partial derivatives of h a (f) 
with respe ct to the f our angles 6s, <fis, §L, 4>l, using 
Eqs. ( 3.15| ) and ( 3.32| ). The partial derivatives of h a (f) 
with respect to the other 5 parameters are given by 
Eqs. (3.14) and (3.19) of @. 



(5.8) 



B. Results 

Using Eq. (4.3) as our model for the signal, we have 
computed the variance-covariance matrix Y? 3 for a wide 
range of masses and angles. We note that since the signal- 
to-noise is high for SMB H mergers, the Fischer matrix 
approximation, Eq. ( |2.8| ), is expected be quite accurate. 

Clearly there is a very large, non-trivial parameter 
space to explore: M, fi, /3, 9s, 4>s, 9~l, 4>l- (The Fisher 
matrix is independent of </> c and t c , and lnD^ just af- 
fects the overall scaling.) Here we will look at just a few 
representative cases. Our results are shown in Table II. 
In all cases we take (3 = 0, (that is, the true value of (3 
is zero, but the best-fit value can be non-zero), and we 
take the initial position and orientation of the detector 
to be 0o = 0, a = 0. We take as our fiducial source a 
binary at z = 1 in a low-density (fio = 0) universe with 
Hq = 75 km/s-Mpc; consequently our fiducial distance is 
D L = Hq 1 = 12.253 x 10 9 yr. The masses listed in Table 
II are the 'true' masses, as they appear in Eq. (5.2), not 
their redshifted versions. All results are for one year of 
data. 

As stated earlier, for M. > 3 x 10 6 M Q , most of the 
signal-to-noise accumulates at the very end of the inspi- 
ral, and therefore the total S/N we calculate with our 
model waveform depends rather sensitively on how one 
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assumes the signal is 'cut-off' as the two bodies merge. 
Therefore the S/N results listed in Table II should not 
be regarded as accurate to better than a factor of ~ 2. 
However we find that our results for Af2s do not depend 
sensitively on this cut-off (for a given source at fixed dis- 
tance). We also find that, for the mass range we looked 
at (10 4 to 1O 7 M0), increasing the observation to include 
longer than the final year did not significantly increase 
the angular resolution. 

The following points emerge from Table II. Unlike the 
case for monochromatic sources, having essentially two 
detectors does substantially increase LISA's angular res- 
olution for SMBH mergers. The angular resolution Af^s 
achievable by detectors I and II combined is roughly 10~ 4 
steradians. The angular resolution depends strongly on 
the masses and the particular angles involved, however. 
Afl$ is roughly in the range 10~ 5 — 10~ 3 steradians for 
1O 5 M < M(l + z) < 10 7 M©. AQ S is generally larger 
for the lower mass black holes (M « 10 4 — 1O 5 M ) 
because the signal-to-noise is generally smaller in these 
cases. 

The angular resolution achievable by detectors I and II 
combined is roughly an order of magnitude better than 
that achievable with detector I alone. Notice this is quite 
different from the case of monochromatic sources, where 
the improvement was only a factor of ~ 2. It seems clear 
that this difference arises because, in the SMBH case, 
the time-scale over which most of the signal-to-noise is 
accumulated is rather shorter than a year. During this 
'effective' integration time the orientation of detector I 
does not change dramatically, and so detector I is effec- 
tively sensitive to only a single polarization. 

The distance determination accuracy ADl/Dl for 
SMBH mergers is roughly in the range 0.1% — 30%, 
with ~ 1% being typical. This is much worse than the 
naive guess of ADl/Dl ~ (S/N)" 1 . Clearly the 'extra 
error' is due to correlations between Dl and the vari- 
ous angles describing the source. The quantity Af2^ = 



2tt 



AJJl l A0 L -(A/2 L A</> L ) 



represents LISA's accuracy 



in determining the orbital plane of the binary. We find 
that large values of ADl/Dl have a strong positive cor- 
relation with large values of Afi^, as one might expect. 

Finally we note that a recent paper by Schilling ]23| ] 
has looked more carefully LISA'a transfer function, and 
concludes that the high-frequency part of the instrumen- 
tal noise is somewhat lower than previously estimated, 
so that the term (3(f) in Eqs. (3.34)~(3.35) should be re- 
duced by some factor < 1.5. Re-running the code with 
this change, we find, makes a negligible correction to the 
values of A fig listed in Table II. 



VI. CONCLUSIONS AND FUTURE WORK 

We have seen that, for SMBH mergers, LISA should 
achieve an angular resolution of (very roughly) ~ 0.3 



square degrees. What are the implications of this re- 
sult for the idea of using such mergers to determine the 
cosmological parameters H , £1 , and A ? Since LISA 
can determine the luminosity distance to the source but 
not its redshift, one clearly needs the reshift of the host 
galaxy or galactic cluster to do cosmology. Since one 
square degree contains ~ 10 4 galaxies, the LISA mea- 
surement alone is clearly insufficient to identify the host 
galaxy or cluster. On the other hand: because the signal- 
to-noise is so large, one will know that a merger is oc- 
curring weeks before the final burst. We have checked 
that, more than a day before the final burst, LISA will 
have achieved most of the angular resolution indicated in 
Table II. Thus one will know very accurately when the 
final burst will occur, and will know to within a degree 
where it will occur. At the right time, every available 
telescope can be trained at the right spot in the sky (as 
happened with the impact of the comet Shoemaker- Levy 
on Jupiter), and, if one is lucky, the merging binary could 
'send up a flare' electromagnetically (t]]. Of course, a flare 
is possible only if there is normal matter involved in the 
collison; e.g., if at least one of the black holes has pre- 
served an accretion disk up to the point of merger. It is 
possible there could be some electromagnetic signal even 
before the final burst. Clearly, these possibilities deserves 
more investigation. 

Finally, we list some ways in which our analysis could 
be improved. First, it would be useful to have a better 
understanding of how the instrumental noise in the dif- 
ferent arms will be correlated; our assumption of 'total 
symmetry' between the different arms was a crude esti- 
mate intended just to get us started. Of course it would 
also be useful to have a better estimate of the confusion 
noise levels, but it may take a working LISAto provide 
that! Second, our treatment of monochromatic sources 
could be improved by doing a full Monte Carlo estimation 
of the errors. The S/N for any monochromatic source is 
likely to be low, it is not clear how reliable our Fischer 
matrix approxiamtion is in this case. Thirdly, one could 
clearly do a much more systematic analysis of the pa- 
rameter space: the parameters in Table II were chosen 
more or less at random. Finally, a better treatment of 
the SMBH case would allow for eccentric orbits and spin- 
induced precession of the binary's orbital plane. Work to 
incorporate these last two effects in our analysis is now 
in progress p|Jl"o[|. 
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TABLE I. LISA's angular resolution Afls (in steradians) 
for monochromatic sources. Results are for a 1-year obser- 
vation, with source strength normalized so that total (i.e., 
detectors I and II combined) S/N = 10. Si/N and AQ s ,i 
are the signal-to-noise and angular resolution achievable by 
detector I alone, for the same source strength. LISA's initial 
position and orientiation are given by ao = <f>o = 0. 
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TABLE II. LISA'S angular resolution for SMBH mergers. All mergers are taken to occur at redshift 2 = 1 and luminosity 
distance Dl = H^ 1 , where we take Ho = 75 km/s-Mpc. Results marked with subscript "I" are for detector I alone; results 
without a subscript represent the signal-to-noise and acuracies achievable using both detectors I and II. LISA's initial position 
and orientiation are given by «o — 4>o — 0. 
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